The Altitude of an Infrared Bright Cloud Feature on Neptune 
from Near-Infrared Spectroscopy 1 
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ABSTRACT 

We present 2.03-2.30 /mi near-infrared spectroscopy of Neptune taken 1999 June 2 (UT) with 
the W.M. Keck Observatory's near-infrared spectrometer (NIRSPEC) during the commissioning 
of the instrument. The spectrum is dominated by a bright cloud feature, possibly a storm or 
upwelling, in the southern hemisphere at approximately 50? S latitude. The spectrum also includes 
light from a dimmer northern feature at approximately 30?N latitude. We compare our spectra 
(A/AA ^2000) of these two features with a simple model of Neptune's atmosphere. Given our 
model assumption that the clouds are flat reflecting layers, we find that the top of the bright 
southern cloud feature sat at a pressure level of 0.14 i !o3 bar, and thus this cloud did not extend 
into the stratosphere (P<~0.1 bar). A similar analysis of the dimmer northern feature gives a 
cloud-top pressure of 0.084±0.026 bar. This suggests that the features we observed efficiently 
transport methane to the base of the stratosphere, but do not directly transport methane to 
the upper stratosphere (P< 10 -2 — 10~ 3 bar) where photolysis occurs. Our observations do not 
constrain how far these clouds penetrate down into the troposphere. We find that our model fits 
to the data restrict the fraction of H2 in ortho/para thermodynamic equilibrium to greater than 
0.8. 

Subject headings: infrared: solar system — planets and satellites: Neptune 
1. Introduction 
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The first hint of Neptune's atmospheric com- 
plexity and variability came when Joyce et al. 
(1977) observed significant changes in Neptune's 
brightness at 1 — 4 /im over the course of ap- 
proximately an Earth year. Pilcher (1977) in- 
terpreted this as the formation and slow dissi- 
pation of an extensive high-altitude cloud. The 
1989 flyby of the Voyager II spacecraft revealed a 
host of time-varying atmospheric features (Smith 
et al. 1989). Even before the Voyager II flyby, 
the development of the Charge Coupled Device 
(CCD) allowed imaging of Neptune at wavelengths 
up to ~1 /im. Several observers looking in the 
0.62 and 0.89 /im methane absorption bands reg- 
ularly found mid-latitude features that were ex- 
tremely bright relative to Neptune's disk (Smith 
1984, 1985; Hammel & Buie 1987; Hammel 1989; 



Hammel et al. 1989; Hammcl 1990). These fea- 
tures are presumably clouds and may be storms or 
large upwellings of material from the troposphere. 
When present, the reflected sunlight from these 
features dominates images of Neptune at methane- 
absorbing wavelengths between 0.6 and 2.5 /zm, 
as shown by many observers. Hubble Space Tele- 
scope (HST) regularly observed such features at 
wavelengths less than 1 fim (Sromovsky et al. 
1995; Hammel et al. 1995; Hammcl & Lockwood 
1997), while ground based observers using conven- 
tional infrared techniques have seen these features 
at 1 to 2.5 /zm (Sromovsky et al. 2001a, b,c). More 
recently, high resolution techniques such as speckle 
imaging and adaptive optics (AO) have been used 
to observe Neptune and these bright features at 
1-2.5 /im (Roddicr et al. 1997, 1998; Roe et al. 
2000; Gibbard et al. 2000; Max 2000). 

Speculation about the nature and origin of 
these phenomena has primarily focused on the 
idea of large upwellings punching through the 
tropopausc resulting in a high column density of 
condensed methane particles. Thus, these fea- 
tures could in part be responsible for transporting 
methane through the cold-trap of the tropopausc 
and loading the stratosphere with methane gas, 
where it is then photolyzed and converted to a va- 
riety of heavier hydrocarbons, eventually forming 
hazes (Baines et al. 1995b; Romani et al. 1993; 
Moses et al. 1995). It is crucial for our under- 
standing of the dynamics and chemistry of Nep- 
tune's atmosphere to know the altitude range to 
which these cloud features reach. 

Hammel et al. (1989) estimated from their 
CCD photometry that the bright features they ob- 
served were due to increases in the number den- 
sity of high stratospheric haze particles. Rod- 
dier et al. (1998) observed Neptune with adap- 
tive optics techniques. They used two narrowband 
filters centered on 1.56 and 1.72 /im, such that 
one filter is centered on a strong methane absorp- 
tion feature while the other filter is outside the 
strong methane absorption. These authors esti- 
mated that the bright features are located near the 
tropopause at pressures on the order of 0.1 bar or, 
possibly, at even higher altitudes. More recently 
Sromovsky et al. (2001c), using IRTF photom- 
etry, found the altitudes of a number of discrete 
cloud features to be between 0.060 and 0.230 bar. 
In this paper we present spectra of two of these 



cloud features. Through comparison with a sim- 
ple radiative transfer model, we used these spec- 
tra to determine precisely the altitude of the cloud 
features. Our best-fit model places the top of the 
bright southern cloud feature that we observed at 
a pressure level of 0.14 bar within an uncertainty 
range of 0.11 to 0.19 bar, while the best-fit for the 
dimmer northern feature puts it at 0.084±0.026 
bar. 

2. Observations and Data Reduction 

We observed Neptune on 1999 June 2 (UT) us- 
ing NIRSPEC, the W.M. Keck Observatory's new 
near-infrared spectrometer, on the Keck II tele- 
scope during the commissioning of the NIRSPEC 
instrument (McLean et al. 1998). This spectrom- 
eter operates over a wavelength range of 0.95- 
5.5 fim in either a low-resolution (R ~ 2000) mode 
or a cross-dispersed high-resolution (R ~ 25,000) 
mode. NIRSPEC is equipped with a 1024 x 1024 
InSb ALADDIN array for spectroscopy, and also 
a slit-viewing camera (SCAM) containing a 256 
x 256 HgCdTc PICNIC array with 07 18 pixels. 
The data presented here are from a single low- 
resolution setting using the NIRSPEC-6 block- 
ing filter, and they cover roughly 2.03-2.30 /jm. 
In low-resolution mode the pixel size in the spa- 
tial direction of the ALADDIN spectral array is 
0.144" /pixel. 

We acquired a series of slit-viewing camera 
(SCAM) images both before and during our spec- 
tral exposures, giving us images of Neptune with 
the spectrometer's slit offset from the disk of Nep- 
tune and overlapping the disk of Neptune. SCAM 
images were taken in pairs, and after the first ex- 
posure of each pair the pointing of the telescope 
was offset by 10" . In the current work we did 
not attempt precise photometry, and therefore our 
processing of the SCAM images is simplistic: We 
subtracted one image of each pair from the other 
for background and bias subtraction. We then 
shifted and coadded the images using the Gaussian 
centroids of Triton and an unidentified star for off- 
set determination. From Triton (apparent diam- 
eter 07126) and the unidentified star the FWHM 
of the point spread function was 0743±0707. 

Figure la shows a SCAM image of Neptune 
taken simultaneously with the spectra presented 
in this work where the disk is bisected by the slit. 



Figure lb shows an unobstructed image taken 25 
minutes earlier. The SCAM image shown in Fig. 
la was taken with the NIRSPEC-6 filter, (1.56- 
2.32 /zm), while the unobstructed SCAM image 
in Fig. lb was taken with the NIRSPEC-7 filter 
(1.84-2.63 /xm). We did not take an unobstructed 
SCAM image in the NIRSPEC-6 filter, and there- 
fore we present the NIRSPEC-7 filter image to 
show Neptune unobstructed by the spectrometer 
slit. Neptune's apparent diameter (at 1 bar level) 
was 2'.'30, the Earth's planetographic sub-latitude 
on Neptune was -28?08, and the solar phase an- 
gle was 1?54. 8 Figure lc shows the orientation 
and scale of Neptune on the images in Fig. la- 
fa. Neptune's brightness along the slit is shown 
in Fig. Id. Comparison of Neptune's brightness 
as a function of position on the slit with that 
of HD201941 (Fig. Id) shows that the projected 
size of the storm is only marginally resolved. The 
local minimum between the features on Neptune 
indicates that we can extract spectra of the two 
separate features relatively cleanly without much 
cross-contamination. 

The spectra presented here come from two 60- 
second exposures taken with a 42'.' x 0'.'380 slit 
starting at 14:06 (UT) on 1999 June 2. The slit 
was aligned parallel with Neptune's north-south 
axis and centered on the bright feature in the 
southern hemisphere. This feature was by far the 
brightest that we observed on Neptune on 1999 
June 2 (UT). The slit also captured light from a 
dimmer feature in Neptune's northern hemisphere. 
Between the two exposures the pointing of the 
telescope was moved ~10" along the direction of 
the slit, so that Neptune fit easily on the slit for 
both exposures. In order to correct for Earth's at- 
mospheric absorption, we observed an A2 spectral 
type star (HD201941) in two 10 second exposures, 
with an offset of the telescope pointing between 
exposures in order to move the star along the slit. 

The reduction sequence consisted simply of sub- 
traction of one star spectral frame from the other 
for bias and background subtraction. In these im- 
ages the spatial and spectral coordinates are dis- 
torted with respect to the rows and columns of 
the detector array. The OH sky emission lines 
in unsubtractcd frames trace out lines of constant 




Fig. 1. — (a) SCAM image of Neptune with spec- 
trometer slit overlying Neptune. This image was taken 
in the NIRSPEC-6 filter (1.558-2.31 /im) at 14:06 UT 
during the first Neptune spectral exposure, (b) SCAM 
image of Neptune unobstructed by slit. This image 
was taken in the NIRSPEC-7 filter (1.839-2.630 /urn) 
at 13:45 UT. No unobstructed images of Neptune were 
taken in the NIRSPEC-6 filter, (c) Schematic show- 
ing the scale and orientation of Neptune in (a) and 
(b). North is up and Neptune appears as it would 
look on the sky. The scale of the schematic is shown 
at lower- left and is the same as in (a) and (b). The 
vertical lines show the approximate location of the slit 
in (a), (d) Brightness of Neptune along the slit, av- 
eraged over all the wavelengths shown in the spectral 
image in Fig. 2b. For comparison to atmospheric see- 
ing, the profile of the star HD201941, averaged over 
the same wavelength range, is also shown. The profile 
of HD201941 gives a seeing FWHM of - 0'.'5 - 0'.'6. 
The FWHM of the unidentified star in the the scam 
images taken simultaneously with the Neptune spec- 
tra was 0'.'4. The narrowness of the FWHM compared 
to the separation of Neptune's northern and southern 
features, along with the clear local minimum between 
the features, leads us to conclude that we have sepa- 
rated the light from these two features relatively well. 



3 See the NASA/JPL Horizons ephemeris program at 
http://ssd.jpl.nasa.gov/horizons.html. 



wavelength and are well fit by straight lines 
Y(x)=A (X)+A 1 (X)x, 



(1) 



where A and Ai are functions of wavelength. 
Meanwhile, the arc of a stellar spectrum traces 
out a line of constant position along the slit and 
is well fit by the function 

X(y) = B Q ( S )+B 1 ( S )y + B 2 (s)y 2 , (2) 

where Bo, Bi, and B2 are all functions of slit po- 
sition (s). The (x,y) position on the array for a 
given slit-position and wavelength (s,A) can then 
be found by: interpolating Ao and Ai for A from 
the numerous OH lines that we fit, interpolating 
Bo, Bi, and B2 from the several stellar spectra fit, 
and finally finding the (x,y) intersection of Eqns. 
1 and 2. By doing this for a grid of wavelengths 
and positions along the slit the data in a spectral 
image is interpolated from (x,y) to (s,A). The final 
step in this rectification process is to apply a Jaco- 
bian correction for the geometric distortion. We 
extracted the stellar spectrum from the rectified 
spectral image using an optimal weighted extrac- 
tion technique that includes a median filter rejec- 
tion algorithm to remove the effects of bad pixels 
and cosmic ray hits. Finally, we divided the ex- 
tracted stellar spectrum by that of Vega (spectral 
type AOV) (Colina, Bohlin, & Castclli 1996b) to 
produce an estimate of the combined atmospheric 
and instrumental transfer function, shown in Fig. 
2a. 

We processed the spectrum of Neptune in a 
manner similar to that applied to the calibration 
star; however, after rectification, but before ex- 
traction, we inserted the additional step of divid- 
ing by the transfer function determined from the 
star. The final rectified spectral image of Neptune 
is shown in Fig. 2b. We extracted the spectra of 
the cloud features in a similar manner as for the 
stellar spectrum, except that we limited the ex- 
tractions to 0'.'6 along the slit centered on each 
feature. The spectrum of the dimmer northern 
cloud feature is shown in Fig. 2c, and the spec- 
trum of the brighter southern feature is shown in 
Fig. 2d. Shown in Fig. 2c and 2d are averages 
from the two separate exposures. Having two sep- 
arate exposures provides a check on our precision. 
The two spectra of the northern feature extracted 
from the two exposures appeared identical except 



for random noise. Similarly, the two southern fea- 
ture spectra were also very nearly identical. 

Navigation on the images of Neptune is difficult 
because light from the cloud so dominates over all 
other features, however the presence of Triton in 
the slit-viewing camera images makes this prob- 
lem significantly easier. By centroiding a gaussian 
on Triton and using the offset from Triton to Nep- 
tune's center given by JPL's Horizons ephemeris, 
we find the center of Neptune. Combining this 
with centroiding a Gaussian on the cloud feature 
we estimate the cloud to be located 0'.'65±0'.'25 
from the center of Neptune's disk at a Neptune 
latitude of -48?±6?. Thus, the cloud lay at a view- 
ing angle 9 of 34? ±17?, where 8 is the angle be- 
tween the normal on Neptune's 'surface' and our 
line of sight. By a similar procedure we estimate 
the observed northern feature to be at a latitude 
of 30? ±13? and a viewing angle of 55?+ 9 L 

3. Atmospheric Model 

Our aim in the work presented here is to mea- 
sure the altitude or pressure level at the top of 
an infrared-bright cloud on Neptune. Towards 
this end we have taken spectra, presented in the 
previous section, over a wavelength range where 
the opacity in Neptune's atmosphere varies signif- 
icantly as a function of wavelength due to H2 col- 
lision induced absorption (H2-CIA) and methane 
absorption. In our model we calculate the pre- 
dicted spectrum as a function of the altitude of 
the top of the cloud and several other param- 
eters described below. We judge the goodness- 
of-fit for each model spectrum using the metric 
XX-fohs(A) - A ■ ImodeiW 2 , where A is chosen in 
each case to minimize the overall sum. The in- 
troduction of the factor A is necessary due to the 
lack of an absolute flux calibration for the observed 
spectrum, I h s (A). 

Our model atmosphere consists of 120 layers 
evenly spaced in Logio(P& a r) from 5.0 to 10 -4 bar. 
We interpolate the temperature and pressure for 
each layer from Lindal (1992). The free param- 
eters in our model arc: the mole fraction of he- 
lium, F# e ; the mole fraction of methane in the 
stratosphere, Fcff 4 , s ; the mole fraction of methane 
in the troposphere, FcH^.t] the fraction of H2 
in ortho/para thermodynamic equilibrium, F eq ; 
the viewing angle, O; and the pressure altitude 



of the top of the cloud in bar, Pf, ar . Around 
the tropopause the fractional methane abundance 
follows the saturation vapor curve, so that the 
methane abundance is never super-saturated. 
Wavelengths of 2.03-2.30 /zm do not probe sig- 
nificantly into the troposphere, and therefore our 
model fit is insensitive to changes in FcH 4 .t- In 
each layer a fraction of the H2, F eq , is distributed 
between ortho and para states according to ther- 
modynamic equilibrium, with the remaining H2 
distributed according to an ortho:para ratio of 
3:1. We calculated the model predicted spectrum 
for each point on a grid of these free parameters. 
The grid points for each parameter are listed in 
Table 1. 

To model collision-induced absorption by hy- 
drogen (H2-CIA) for H2— H2 and H2— He collisions 
we use the FORTRAN routines of A. Borysow 
(Borysow 1991, 1993; Zheng & Borysow 1995; Bo- 
rysow et al. 1989a,b). 9 Although both 0-1 and 0-2 
transitions are included in our model, for wave- 
lengths of 2.1-2.3 /im only the 0-1 transition is 
relevant. 

Accurate modeling of methane absorption 
across the near-infrared spectrum is extremely 
difficult due to the enormous number of individ- 
ual lines and huge variation in line strength. We 
apply the correlated k distribution method as de- 
scribed in Lacis & Oinas (1991) and on p. 230 of 
Goody & Yung (1989). We use the H2-broadened 
methane ^-coefficients of Irwin et al. (1996) since 
Neptune's atmosphere is primarily H2. These co- 
efficients are for 5 cm -1 wide bins and this places 
a limit on the spectral resolution of the model. 

We ignore all scattering processes, except re- 
flection from the top of the cloud. Although sig- 
nificant at shorter wavelengths, Rayleigh scatter- 
ing is negligible at wavelengths of 2.0 to 2.3 /im. 
Light reflected from the top of the cloud domi- 
nates all other sources which might contribute to 
our cloud spectrum, such as scattered light from 
stratospheric hydrocarbon hazes. Therefore, the 
only scattering process that we include is reflec- 
tion from the top of the cloud, which we model as 
a flat reflecting layer. The reflectivity of the top 
of the cloud, or alternatively, the combined optical 
depth, scattering phase function, and single scat- 
tering albedo, are irrelevant given that we fit the 



shape of the spectrum, not the absolute flux level. 
Further, we assume that the optical depth, scat- 
tering phase function, and single scattering albedo 
do not vary significantly over the wavelength range 
of our spectrum (2.03 to 2.3 /im). Using the so- 
lar spectrum of Colina, Bohlin, & Castelli (1996a), 
the model produces the estimated spectrum of the 
cloud as a function of the pressure level of the top 
of the cloud. 

4. Results of Model Fit 

The spectral resolution of the model described 
in the previous section is limited by the k- 
cocfficicnts of Irwin et al. (1996) and is coarser 
than the spectral resolution of our observed spec- 
tra. Therefore, we bin the observed spectra in 
wavelength to achieve a resolution as nearly iden- 
tical as possible to the model spectra. We restrict 
the model fitting to the wavelength range 2.08- 
2.25 (im in order to avoid a large telluric CO2 
band at < 2.08/xm and a series of sharp methane 
features at > 2.25/im that are poorly represented 
by the methane coefficients in the model. 

For the bright southern feature at a viewing an- 
gle of = 34?, we obtain our best-fit for a cloud 
top at 0.14 bar, F C h 4 ,s = 0.0017, F He = 0.22, 
and F eq = 1.0. For the dimmer northern feature 
at a viewing angle of 9 = 55?, the best fit is for a 
cloud-top pressure of 0.084 bar, Fch 4 ,s = 0.0017, 
F He = 0.22, and F eq = 1.0. These best-fit model 
spectra are superposed on the observed spectra in 
spectrum (1) of Fig. 3a and spectrum (1) of Fig. 
3b. The two parameters that we can best con- 
strain are Pbar at the top of cloud and the frac- 
tional equilibrium of H2 in ortho /para equilibrium, 
F eq . The model spectra do not fit the observations 
for cloud-top pressures outside the range of 0.11- 
0.19 bar for the bright southern feature (see spec- 
tra (2) and (3) of Fig. 3a), nor outside the range 
of 0.058-0.110 bar for the dimmer northern feature 
(see spectra (2) and (3) of Fig. 3b). As shown in 
Fig. 4, reasonable model fits to both the northern 
and southern spectra require F eq >~0.8. We find 
that our data do not constrain significantly F# e 
and F C h 4 ,s- 

Due to the low spatial resolution of our data 
there is significant uncertainty in the viewing an- 
gle for both features, = 34? ± 17? for the bright 

+9° 

3 Available at: http://www.astro.ku.dk/~aborysow/programs/index.M4illthcrn feature and = 55? A for the dim- 



Table 1 
Model parameter grid points. 



Parameter 



Values used in model 



Fife 0.08, 0.10, 0.12, 0.14, 0.16, 0.18, 0.20, 0.22 

? CH 4 ,a b 2. x 10~ 5 , 8. x 10~ 5 , 1.8 x 10~ 4 , 3.5 x 10 -4 , 7.0 x 10~ 4 , 1.05 x 10~ 3 , 1.7 x 10" 

7 CH 4 ,',t c 0.022 

F e , d 0.0, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 

6 17?, 34?, 51? 



120 layers evenly spaced in Logio(P) from 5.0 to 10 bar 



"■Conrath et cd. (1993) found F He = 0.15 

Baincs & Hammcl (1994) found Fch 4)S — 3.5 X 10~ with a maximum uncertainty range of 
2.5 x 10~ 5 to 1.7 x 10~ 3 . 

c Baincs et al. (1995a) found ^CH^.t — 0.022. The fit of our model to data is not sensitive to 
FcH 4 ,t- 

Baincs & Smith (1990) found F eq — 1.0, with a minimum allowed value of 0.85. 





Fig. 2. — (a) Observed spectrum of A2 star HD201941 
divided by the reference spectrum of AOV star Vega, 
providing an estimate of the atmospheric and instru- 
mental transmission functions, (b) Rectified spectral 
image of Neptune. The spectral region of 2.03 to 2.05 
^m is omitted in order to show an intensity-bar and 
the spatial size of Neptune's diameter along the slit, 
(c) Extracted spectrum of the dimmer northern fea- 
ture, d) Extracted spectrum of the brighter southern 
feature. Note the short vertical tick marks that show 
methane absorption lines identified from the labora- 
tory spectrum of McKellar (1989). The vertical scale 
in (c) is exaggerated by a factor of 11 with respect to 
the vertical scale in (d). 



mer northern feature. Decreasing O to 17? for the 
bright southern feature pushes the best-fit cloud- 
top pressure to 0.16 bar, while increasing to 51? 
changes the best-fit cloud-top pressure to 0.12 bar. 
Similarly for the northern feature, decreasing O 
to 51? moves the best- fit cloud-top pressure to 
0.092 bar, while increasing O to 64? shifts the best- 
fit cloud-top pressure to 0.076 bar. In all these 
cases the best- fit parameters include F eq = 1.0 
and F CHitS = 0.0017. 

5. Errors and Uncertainties 

At this point it is worthwhile to make a brief 
discussion of how the errors and uncertainties in 
our observations and model fitting could affect our 
results with respect to cloud-top pressure and F eq . 

On the observing side, we are much more con- 
cerned with systematic errors, for instance artifi- 
cial slopes across the entire spectrum, than with 
random errors in the spectra of HD201941 and 
Neptune. Since we are fitting the model to 73 
wavelength bins, random errors from bin to bin 
will tend to cancel out and not bias the model fit. 
There are several possible sources of systematic 
errors on the observing side; the three of great- 
est concern relate to alignment on the slit and the 
method of atmospheric correction. Misalignment 
of the slit on the star would redden the spectrum 
and possibly introduce a bias in the final model 
fitting. This is less of an issue on an extended 
source such as the clouds on Neptune. By looking 






Fig. 3. — Constraint on the cloud top pressure, (a) 
In each case the solid dark line is the binned observed 
spectrum of the bright southern feature. The overplot- 
ted grey lines are the model spectra: (1) The over- 
all best fit model is for P ftar =0.144 bar, Fch 4 , s = 
1.7 x 1(T 3 , F He =0.22, F e ,=1.0, and 6 = 34?. (2) 
When the cloud top is raised to 0.092 bar, the best fit 
is poor and requires Fch 4 ,s = 1.7 x 10~ 3 , Fhb = 0.08, 
F eq = 1.0, and 9 = 34?. (3) Similarly, when the 
cloud top is lowered to 0.227 bar, the best fit is 
poor. In this case F C h 4 ,s = 2. x 10 -5 , F He — 0.22, 
F eq = 1.0, and = 34?. (b)In each case the solid 
dark line is the binned observed spectrum of the dim- 
mer northern feature. The overplotted grey lines are 
the model spectra: (1) The overall best fit model is 
for P b ar=0.084 bar, F C h 4 , s = 1.7 x 10 -3 , Fjj e =0.22, 
F e ,=1.0, and = 55?. (2) When the cloud top is 
raised to 0.053 bar, the best fit is poor and requires 
Fch 4 , s = 1.7 x 10" 3 , F He = 0.08, F e , = 1.0, and 
= 55?. (3) Similarly, when the cloud top is low- 
ered to 0.132 bar, the best fit is poor. In this case 
F CHi ,s = 1.05 x 10" 3 , F He = 0.22, F e , = 0.9, and 
= 55? . (c) The temperature-pressure profile of Lin- 
dal (1992). Also plotted are the cloud top pressures of 
the model spectra in Fig. 3a. The short vertical line 
on the far left shows the narrow range of pressures to 
which our observations restrict the top of the cloud. 



at multiple stellar spectra we estimate that this 
source of error introduces at most a one to two 
percent slope from 2.08 to 2.25 /im. 

In applying the atmospheric and instrumental 
transmission correction with HD201941 there are 
two more potential sources of systematic error. 
The first is that to find the atmospheric trans- 
mission function we divided HD201941 by a spec- 
trum of Vega, and the second is that HD201941 
was not observed at exactly the same airmass and 
time as Neptune. While Vega is an A0V star, 
HD201941 is listed as an A2 star in the SIMBAD 
database. 10 In order to estimate the maximum 
slope bias that this stellar mis-match could in- 
troduce, we compared blackbody curves. Drilling 
& Landolt (2000) give the T e// for an A0V star 
as 9790 K and for an A2V star as 9000 K. This 
difference suggests a slope error of 0.29 percent 
from 2.08 to 2.25 pm. The spectra of HD201941 
were taken at an airmass 1.05, while the Neptune 
spectra were taken at airmass 1.28. In order to 
minimize the influence of this on the model fitting 
we excluded wavelengths shortward of 2.08 /im to 
avoid a large CO band. To investigate what bi- 
ases and slopes this mismatch in airmass could 
introduce we used the ATRAN (Lord 1992) model 
atmospheric transparency spectra available on the 
Gemini Observatory website. 11 While we did not 
examine transparency spectra for our exact air- 
masses, the slope difference introduced by observ- 
ing at airmass 1.0 versus 1.5 across our spectral 
range of interest would be 1.6 percent. 

Each of these possible slope errors discussed 
above is less than 2 percent. To show that even a 
fortuitous addition of all these slope errors in one 
direction would not change our primary results we 
artificially introduced slope errors of ±10 percent 
to our final observed spectra and refit the model. 
This had no effect on results concerning F eq and 
at most shifted the best-fit pressure level of the 
cloud top by one level in our model, to 0.12 bar 
in the —10% case for the bright southern feature 
and to 0.09 bar in the +10% case for the dimmer 
northern feature. 

While there are numerous small ways in which 
the model may be inaccurate, for instance if the 
temperature-pressure curve is not exactly correct 



10 Available at http://simbad.u-strasbg.fr. 

11 Sec http://www.gemini.edu/sciops/telcscope/telIndex.html. 
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Fig. 4. — Constraint on the fraction of H2 in or- 
tho/para equilibrium (F eq ). (a) In each case the solid 
dark line is the binned observed spectrum of the bright 
southern feature. Holding F eq fixed at 0.6, 0.8, 1.0, 
the best fits of model to data are shown. In all 
three cases the best fit included P(, ar =0.144 bar and 
Fcfl4, s 1.7xl0 -3 . Incases (I) and (3) F/j e =0.22, while 
in case (2) F_H e =0.20, which is not a significant differ- 
ence, (b) In each case the solid dark line is the binned 
observed spectrum of the dimmer northern feature. 
Holding F eq fixed at 0.6, 0.8, 1.0, the best fits of model 
to data are shown. In all three cases the best fit was 
Pi, ar =0.144 bar, F C // 4 , S 1.7 x HT 3 , and F ffe =0.22. 



for the location of the cloud features, the two ma- 
jor sources of uncertainty in the model are the 
methane ^-coefficients of Irwin et al. (1996) and 
the assumption of a flat reflecting cloud layer. For 
Neptune's atmosphere we are forced to extrapolate 
the methane ^-coefficients to much colder temper- 
atures than the temperatures of the laboratory 
measurements on which they are based. While 
the accuracy or inaccuracy of this extrapolation 
is difficult to judge, the independence of best-fit 
cloud-top pressure and F eq from methane concen- 
tration Fch 4 ,s gives us confidence in our results. 
For ease we assume in our model that the cloud 
top is a flat reflecting layer, however due to particle 
scattering properties the reflectivity of the cloud 
may vary with wavelength and the 'top' of the 
cloud is almost certainly somewhat extended. Our 
best-fit model for the bright southern feature (see 
spectrum 1 of Fig. 3a) is systematically slightly 
off from the observed spectrum. This is most eas- 
ily seen at wavelengths shortward of ~2.16 /im 
where H2 absorption dominates. This is sugges- 
tive that while a flat reflecting layer at 0.144 bar 
does not fit perfectly, a combination of reflectance 
from pressures slightly higher to slightly lower 
than 0.144 bar might result in a better fit, which 
is exactly what one would expect if the cloud-top 
were somewhat extended. In fact, a more detailed 
method of modeling would be to view the com- 
puted model spectra for all the pressure levels in 
the model as a basis set and to construct the best- 
fit model spectrum from a linear combination of 
all the spectra from different levels. However, one 
shortfall to that approach is that it implicitly as- 
sumes single scattering, ignoring multiple scatter- 
ing between layers. A more complete modeling 
approach must include the multiple scattering be- 
tween layers as well, which is beyond the scope 
of the current paper. Variation in reflectivity as 
a function of wavelength may also play a role in 
these slight discrepancies between model and data. 

6. Conclusions 

By comparing a near-infrared spectrum with 
the predictions of a simple transmission model 
we determined the pressure level at the top of 
an infrared-bright tropospheric cloud on Neptune. 
We find a best-fit of model to data for a cloud- 
top pressure level of 0.14 bar within a maximum 



allowed range of 0.11 to 0.19 bar for the bright 
southern feature that we observed. We found the 
dimmer northern feature to sit slightly higher in 
the atmosphere at 0.084 bar within a maximum al- 
lowed range of 0.058 to 0.11 bar. Our work places 
no limit on the pressure at the bottom of the cloud. 
Our results further restrict the fraction of H2 in or- 
tho/para equilibrium to greater than 0.8, and our 
best- fits consistently put this fraction at 1.0. This 
is in agreement with the work of Baines & Smith 
(1990) who found the same results, but from a dif- 
ferent technique, measuring the equivalent widths 
of the 4-0 S(0) and S(l) transitions between 0.6 
and 0.7 /im. Our results do not constrain the frac- 
tional abundance of methane in the stratosphere, 
nor the fractional abundance of helium. 

Our primary result is the tight constraint we 
place on the pressure at the top of the cloud. By 
constraining the cloud-top to pressures around the 
tropopause, we show that the cloud, possibly a 
storm or upwelling, does not extend significantly 
into the stratosphere. If the cloud is made up of 
condensed methane particles brought up from be- 
low, then the mechanism by which this cloud was 
formed appears to be efficient at bringing methane 
to near the top of the troposphere, but, at least at 
the time we observed, the mechanism was not act- 
ing as an efficient method of transporting methane 
to the upper stratospheric levels where ultraviolet 
photolysis occurs (P> 10~ 2 — 10~ 3 bar). 

In the current paper we present a measurement 
of the altitude of two cloud features at a single 
time. We expect longer-term observations of mul- 
tiple infrared-bright features will find that most 
reach only to approximately the tropopause, as 
in the case presented here, but occasional fea- 
tures may reach far into the stratosphere (P> 
0.01 bar) and thus would provide an extremely ef- 
ficient method of transporting methane to the up- 
per stratosphere for photolysis. We are currently 
undertaking such a program of observations using 
NIRSPEC coupled to the Keck Adaptive Optics 
system (Wizinowich et al. 2000) to achieve simul- 
taneous high-spatial and high-spectral resolution. 
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